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Abstract: We propose a first implementation of tlie integrand- reduction metfiod for two- 
<^ \ loop scattering amplitudes. We show that the residues of the amplitudes on multi-particle 

cuts are polynomials in the irreducible scalar products involving the loop momenta, and 
that the reduction of the amplitudes in terms of master integrals can be realized through 
polynomial fitting of the integrand, without any apriori knowledge of the integral basis. We 
discuss how the polynomial shapes of the residues determine the basis of master integrals 
appearing in the final result. We present a four-dimensional constructive algorithm that 
we apply to planar and non-planar contributions to the 4- and 5-point MHV amplitudes 
in A/" = 4 SYM. The technique hereby discussed extends the well-established analogous 
method holding for one- loop amplitudes, and can be considered a preliminary study towards 
the systematic reduction at the integrand-level of two- loop amplitudes in any gauge theory, 
suitable for their automated semianalytic evaluation. 
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1. Introduction 

Scattering amplitudes constitute the core of the perturbative structure of quantum field 
theories. The recent development of novel methods for computing them has been highly 
stimulated by a deeper understanding of the multi-channel factorization properties natu- 
rally emerging as reactions of the amplitudes under deformations of the kinematics in the 
complex plane dictated by on-shell [1,2] and generalized unitarity-cut conditions [3,4]. 

Analyticity and unitarity of scattering amplitudes [5] have then been strengthened by 
the complementary classification of the mathematical structures present in the residues 



at the singularities, better understood after uncovering a quadratic recurrence relation 
for tree-level amplitudes, the so called BCFW-recursion [2], its link to the leading singu- 
larity of one-loop amplitudes [4], and a relation between numerator and denominators of 
one-loop Feynman integrals, yielding the multipole decomposition of Feynman integrands, 
stronghold of the by-now known as OPP method [6]. 

These new insights, which stem from a reinterpretation of tree-level scattering within 
the twistor string theory [7] , have catalyzed the study of novel mathematical frameworks 
in the more super symmetric sectors of quantum field theories, such as dual conformal sym- 
metries [8], grassmanians [9], Wilson-loops/gluon-amplitudes duality [10], color /kinematic 
and gravity/gauge dualities [11,12], as well as on-shell [1,2,13,14] and generalised unitarity- 
based methods [3,4,15,20], and more generally the breakthrough advances in automating 
the evaluation of multi-particle scattering one-loop amplitudes, as demanded by the exper- 
imental programmes at hadron colliders. 

In general, when a direct integration of Feynman integrals is prohibitive, the evalua- 
tion of scattering amplitudes beyond the leading order is addressed in two stages: i) the 
reduction in terms of an integral basis, and ii) the evaluation of the elements of such a 
basis, called master integrals (Mi's). 

At one-loop, the advantage of knowing apriori that the basis of Mi's is formed by 
scalar one-loop functions [16], as well as the availability of their analytic expression [17], 
allowed the community to focus on the development of efficient algorithms for extracting 
the coefficients multiplying each Mi's. Improved tensor decomposition [18], complex in- 
tegration and contour deformation [19], on-shell and generalised unitarity-based methods, 
and integrand-reduction techniques [6, 20-22] led to results which only few years ago were 
considered inconceivable, and to such a high level of automation [21, 23] that different 
scattering processes at the next-to-leading order accuracy can be handled by single, yet 
multipurpose, codes [24-28]. 

At higher-loop, and in particular at two-loop to begin with, the situation is different. 
The basis of Mi's is not known apriori. Mi's are identified at the end of the reduction 
procedure, and afterwards the problem of their evaluation arises. The most used multi- 
loop reduction technique is the well-known Laporta algorithm [29], based on the solution of 
algebraic systems of equations obtained through integration- by-parts identities [30]. The 
recent progress in evaluating amplitudes beyond one-loop has been necessarily accompa- 
nied by the improvement of mathematical methods dedicated to Feynman integrals, such 
as difference [29,31] and differential [32] equations, Mellin-Barnes integration [33], asymp- 
totic expansions [34], sector decomposition [35], complex integration and contour deforma- 
tion [36] - to list few of them. 

In this paper we aim at extending the combined use of unitarity-based methods and 
integrand-reduction, in order to accomplish the semianalytic reduction of two-loop am- 
plitudes to Mi's. The use of unitarity-cuts and complex momenta for on-shell internal 
particles turned unitarity-based methods into very efficient tools for computing scattering 
amplitudes. These methods exploit two general properties of scattering amplitudes, such 
as analyticity and unitarity: the former granting that amplitudes can be reconstructed 



from the knowledge of their (generahsed) singularity-structure; the latter granting that the 
residues at the singular points factorize into products of simpler amplitudes. Unitarity- 
based methods are founded on the underlying representation of scattering amplitudes as a 
linear combination of Mi's, and their principle is the extraction of the coefficients entering 
in such a linear combination by matching the cuts of the amplitudes onto the cuts of each 
MI. 

In the past years, general criteria to determine the Mi's of arbitrary problems have 
been investigated [37], and very recently, a minimal basis for two-loop planar integrals 
has been identified [38], through an ameliorated solution of systems of equations involving 
integration- by-parts identities and supplementary Gram-determinant relations. 

Cutting rules as computational tools have been introduced at two-loop in the context of 
super symmetric amplitudes [39] and later applied to the case of pure QCD amplitudes [40]. 
The use of complex momenta for propagating particles to fulfill the multiple cuts of two- 
loop amplitudes has been proposed for extending the benefits of the one-loop quadruple-cut 
technique, to the octa-cut [41] and the leading singularity techniques [42], as well as the 
method of maximal cuts [43], all indicating the possibility of "reducing the computation of 
multi-loop amplitudes to the computation of residues (which end up being related to tree- 
am,plitudes) and to the solution of linear systems" [42] - an idea we elaborate on hereby. 

The multi-particle pole decomposition for the integrands of arbitrary scattering ampli- 
tudes emerges from the combination of analyticity and unitarity with the idea of a reduction 
under the integral sign. 

The principle of an integrand-reduction method is the underlying multi-particle pole 
expansion for the integrand of any scattering amplitude, or, equivalently, the relation be- 
tween numerator and denominators of the integrand: a representation where the numerator 
of each Feynman integral is expressed as a combination of products of the corresponding 
denominators, with polynomial coefficients. 

The key element in the integrand-decomposition is represented by the shape of the 
residues on the multi-particle pole before integration: each residue is a (multivariate) 
polynomial in the irreducible scalar products (ISP's) formed by the loop momenta and either 
external momenta or polarization vectors constructed out of them; the scalar products 
appearing in the residues are by definition irreducible, namely they cannot be expressed 
in terms of the denominators of the integrand - otherwise mutual simplifications may 
occur and the notion of residues to a specific set of vanishing denominators would become 
meaningless. 

The polynomial structure of the multi-particle residues is a qualitative information that 
turns into a quantitative algorithm for decomposing arbitrary amplitudes in terms of Mi's 
at the integrand level. In the context of an integrand-reduction, any explicit integration 
procedure and/or any matching procedure between cuts of amplitudes and cuts of Mi's is 
replaced by polynomial fitting, which is a simpler operation. 

Decomposing the amplitudes in terms of Mi's amounts to reconstructing the full poly- 
nomiality of the residues, i.e. it amounts to determining all the coefficients of each poly- 
nomial. 



The main goal of this paper is to outhne guiding criteria to constrain the polynomial 
form of the residue on each multiple-cut of an arbitrary two-loop amplitude. Unlike the 
one-loop case, where the residues of the multiple-cut have been systematized for all the 
cuts, in the two-loop case, their form is still unknown. Their existence is a prerequisite for 
establishing a relation between numerator and denominators of any two-loop integrand. 
Their implicit form can be given in terms of unknown coefficients, which are determined 
through polynomial fitting. As in the one-loop case, the full reconstruction of the poly- 
nomial residues is engineered via a projection technique based on the Discrete Fourier 
Transform [44], and requires only the knowledge of the numerator evaluated at explicit 
values of the loop momenta as many times as the number of the unknown coefficients. 

Another feature of the integrand-reduction algorithm we are describing is that the 
determination of the polynomial form of the residues amounts to choose a basis of Mi's, 
which does not necessarily need to be known apriori. In fact, as we will see, each ISP 
appearing in the polynomial residues is the numerator of a potential MI which may appear 
in the final result (other than the scalar integrals). We remark that the set of Mi's which 
will emerge at the end of the integrand-reduction (as well as after applying any unitarity- 
based methods) is not necessarily the minimal set of basic integrals. Integration-by-parts 
identities, Lorentz-invariance identities, as well as Gram-determinant identities may not 
be detected in the framework of a cut-construction, and therefore constitute additional, 
independent relations which can further reduce the number of Mi's which have to be 
evaluated after the reduction stage. 

We define the m— fold cut of a diagram as the set of on-shell conditions corresponding 
to the vanishing of m denominators present in that diagram. As in the maximal-cut 
method [43], we do not cut additional denominators which might arise from cutting one- 
loop sub-diagram. 

We outline the driving principles for the cut-construction of the residues at the m— fold 
cut, and the use of self-consistency checks that ensure the correctness of the reconstructed 
polynomials, namely after the determination of their unknown coefficients. These checks, 
called local and global (N = A^)-tests, are analogous to the tests employed in the one-loop 
integrand-reduction [21,23], and monitor, respectively, the completeness of the polynomial 
residues, and the correctness of the final decomposition formula. 

The values of the loop momenta used for the numerator sampling are chosen among the 
solutions of the corresponding m— fold cut. 

We verify the integrand-reduction algorithm by applying its procedures to planar and 
non-planar contributions to the 4-point MHV [39] and 5-point MHV [47] amplitudes in 
A/" = 4 SYM, and derive an expression for the non-planar pentacross diagram in terms of 
master integrals. 

This work can be considered as a first building block of a new technique, that once 
developed in toto, would allow for the semianalytic reduction of multi-loop amplitudes. 



2. Four-dimensional Reduction Algorithm 

The reduction method hereby presented extends to two-loop the integrand-reduction pro- 
cedures originally elaborated for arbitrary one-loop scattering amplitudes [6,20]. 

An arbitrary two-loop n-point amplitude in the dimensional regularization scheme can 
be written as 



An = f d^-^'q f d^-^'k A{q,k) , 



A{q,k)= ^^^'^^ , (2.1) 

Di = (uiq + j3ik + pif - mf, ai,(3i G {0, 1} 

where e = (4 — d)/2, and d is the continuous-dimensional parameter. Extra-dimensional 
components of the loop momenta can be parametrized as pseudo-mass variables, A^ and 
Afc, one for each loop momenta, according to the following scheme, 

^4-2e = 44 + i\l^ , ql-2e = qI - >^l , (2-2) 

h-2e = h + i^klb , kl_2e = kj - Xl , (2.3) 

'4-2e„ / j4-2er _ /" j-2e \ /" j-2e x , f j4„ /" ,4 



d'^-^'q / d^-^'k = / d-''Xg / d-''Xk / d\4 / C/U4 . (2.4) 

As in the one-loop case, A^ and A^ would appear as additional variables in the polynomial 
residues, and, consequently, would be responsible for the appearance of Mi's in higher 
dimensions [40]. The following discussion is limited to a purely four-dimensional reduction, 
in which the loop variables q^ and k^ are defined in four dimensions. 

2.1 Integrand Decomposition 

The stronghold of the integrand-reduction of an arbitrary two-loop n-point amplitude is 
the decomposition of the numerator N{q, k) in terms of denominators Di for i = 1, . . . , re. 
Following the same pattern as in the one-loop case, a plausible ansatz reads, 

n n n n 

N{q,k)= J2 Ai,_i,iq,k) II Z)ft+ ^ A,,,...,i,((?,/c) J] Dh + 

Ji«ig hj^ii,...,is ii«J7 hj^ii,...,!-/- 

n n 

+ ...+ ^ Ai,^,,{q,k) n Dh, (2.5) 

where ii « is stands for a lexicographic ordering ii < 12 < ■ ■ ■ < ir < is, and where the 
A's are functions depending on the loop momenta. By using the decomposition ( |2.5| ) in 
Eq.(^I^), the multi-pole nature of the integrand of an arbitrary two-loop re-point amplitude 
becomes manifest, 

Af TN Y^ Aj jg(g,fc) Y^ Aij^ i^{q,k) Y^ AjiJ2(g, A;) 

^^''''^= ^ D-D—:d-^ ^ A,A....A. +---+ ^ A,A. • ^'-'^ 



The above expression, upon integration, yields the decomposition of the amphtude in 
terms of Master Integrals (Mi's), respectively associated to diagrams with 8-, 7-, . . . , 2- 
denominators, namely down to the products of two 1-point functions (one tadpole for each 
loop). In Eq.(2.6) each function A{q,k) parametrizes the residue of the amplitude on 



the multi-particle cut that corresponds to the set of vanishing denominators it is sitting 
on. The four-dimensional decomposition in Eqs.( p.5|j2.6| ) begins with 8-denominator terms. 
They correspond to the maximal singularities of two-loop amplitudes in four dimensions, 
accessed by freezing both integration momenta with the simultaneous vanishing of eight 
denominators. We expect that for dimensionally regulated amplitudes, due to the presence 
of additional degrees of freedom, Xg and Afc, an extended integrand-decomposition formula 
should hold, where higher-denominator functions are accommodated. Extension of the 
presented method to two-loop amplitudes in dimensional regularization will be the subject 
of a future work. 

2.2 Residues 

We define the m— fold cut of a diagram as the set of on-shell conditions corresponding to 
the vanishing of m denominators present in that diagram. The calculation of a generic 
scattering amplitude amounts to the problem of extracting the coefficients of multivariate 
polynomials, generated at every step of the multiple-cut analysis. In fact, we will see that 
each A{q, k) is polynomial in the scalar products of the loop momenta with either external 
momenta or polarization vectors constructed out of them. These scalar products cannot 
be expressed in terms of the denominators Dj, and therefore are defined irreducible scalar 
products (ISP's). In the case of the residue to an ?7i— fold cut, Aj^^...^j^, with m < 8, 
the ISP's correspond to the components of the loop momenta not frozen by the on-shell 
conditions; the eightfold-cut conditions of a two-loop amplitude freeze completely both 
integration momenta, like the quadruple-cut in the one-loop case. 

The polynomial form of Aj^^..._j^ depends on the number of independent external mo- 
menta of the n-point diagram identified by Djj . . . Di^ . In an n-point diagram, due to 
momentum conservation, only (n — 1) external momenta are independent. In four dimen- 
sions, there can be at most four independent external momenta. Therefore, despite the 
number of loops, we can trivially cast re— point amplitudes in two groups according to 
whether n is larger than 4 or not: 

• In the case of n— point diagrams with re > 5, four (out of re) external momenta can 
be chosen to form a real four-dimensional vector basis. 

• In the case of re— point diagrams with n < 4, we can choose only up to three (out 
of re) independent external momenta. But they are not sufficient, and additional 
elements, orthogonal to them, have to be taken into account to complete the four di- 
mensional basis. In this case, one can use complex polarization vectors as orthogonal 
complement, to form a complex basis. 

• Each rre- fold cut will be characterized by one of these two kinds of basis, and the 
loop momenta will be decomposed along the vectors forming it. 



• The residue of an m— fold cut, Ajj^,,.^j^, is polynomial in the components of the 
loop momenta, and therefore it is polynomial in the ISP's constructed from the loop 
momenta and the elements of the basis used to decompose them. 

• The polynomial form of the residue Aj^^...^j^ determines the Mi's potentially appear- 
ing in the final decomposition. 

2.2.1 Polynomial Structures and Master Integrals 

Let's consider q and k as the solutions of the m-fold cut identified by the vanishing of 
Di^ . . .Di^. We can decompose the loop momentum q in terms of the basis {tj} and k 
along the basis {cj}, 

g^ = -Po + ^I'^f + ^'i'^2 + ^^STs* + X^T^ , (2.7) 

k" = -r^ + yiel + yse^ + yse^ + y^e^ . (2.8) 

where po and tq are combinations of external momenta. 

For simplicity, let us assume that Aj^^...^j^ is a one-dimensional polynomial. The mul- 
tivariate extension follows the same principles. We can write Aj^^,,,^jj^ in terms of scalar 
products (cj • {k + ro)) as 

A,,_,^=Y,Cj{e^-{k + ro)y . (2.9) 

j 

Since Aj^^...^j^ is the polynomial sitting on the denominators Di^ . . . Di^, the integral ex- 
pression, 

d\ f d^k -^li-^ = ^cJ d\ f d^k (fi_(^i±^ , (2.10) 

J Ui-^ . . . JJi^ J J JJi-^ . . . Ui^ 

can generate Mi's. In general, there are two possibilities which can be encountered: 

A/A%i^^^O; (2^11) 

,4„ / ,4, {ei-{k + rQ)y 



ii) / d\ / d^k -- -'^ -^" = . (2.12) 

In case (i), the natural integral basis associated to the set of denominators Djj . . . Di^ is 
determined by the ISP (cj • (k + ro)). According to the power j, we identify the following 
Mi's: 

i = : f d'^q f d^k ^—-— (scalar MI) (2.13) 

J J ^ii ■ ■ ■ J-Jim 

j = l: f d\ j d'^k (g^-(fe + ^o)) (linear MI) (2.14) 

j = 2: [dS f d^k ^""n " ^^ ^n°^^' (quadratic MI) (2.15) 

J J J^ii ■ ■ ■ J^im 



Option (a) may occur in the case of an ttt,— fold cut of an n-point diagram with n < 4, where 
the vector ef could be orthogonal to the independent vectors of the diagrams, yielding the 
vanishing of the integral in Eq.( 2.12| ). Borrowing the terminology from the one-loop case, 



the ISP is spurious and does not generate any MI, other than the scalar one. 

We observe that in the one-loop case the polynomial residues, from the single-cut up 
to the quadruple-cut, are written in terms of spurious ISP's. Only the pentuple-cut could 
have been written in terms of non-spurious ISP's. However, at the one-loop level there 
is no room for non-spurious ISP, because any scalar product between the loop variables 
and external momenta is reducible. In fact, this is another reason for the one-loop 5-point 
amplitude in four dimensions to admit a complete reduction in terms of 4-point integrals. 

2.2.2 Numerator Sampling 

As in the one-loop integrand-reduction, once the polynomial form of the residues is given 
in terms of unknown coefficients, the numerator decomposition in Eq.(p.5|) becomes the 
engine for determining them. 

Since the function on the l.h.s of Eq.( |2.5D , namely the numerator N{q,k) of an arbitrary 
two-loop diagram is a known quantity, and the r.h.s of Eq.( |2.5D has a canonical represen- 
tation in terms of A-polynomials and denominators Dj^ . . . Di^ , the determination of the 
unknown coefficients sitting in each of the Aj^^...^j. can be achieved by solving a system 
of linear equations. This set of equations is obtained evaluating the l.h.s and the r.h.s of 
Eq.(^]^) for explicit values of the loop variables, as many times as the number of unknown 
coefficients to be determined. The numerator sampling can be performed either numeri- 
cally or analytically. 

Choosing the solutions of the m— fold cuts as sampling values is very convenient, because 
it allows a direct determination of the coefficients appearing in the corresponding residue 
^ii,...,im- Also, the polynomial fitting can be better achieved by a projection technique 
based on the Discrete Fourier Transform [44] , instead of inverting a system. 

2.2.3 Momentum Basis 

Once the polynomial form of each Aj^^..._j^ is established, in terms of spurious or non- 
spurious ISP's, one can use any basis to decompose the loop momenta q and k fulfilling the 
on-shell conditions Di^ = . . . = Di^ = 0. We find it convenient to use always a complex 
four dimensional basis [6], which directly exposes the spurious character of ISP's, should 
they be present. 

For each cut, we decompose the loop momenta q and k, by means of two specific basis 
of four massless vectors. We begin with the decomposition of the loop momentum g'^, 
flowing through an on-shell denominator carrying momentum (q + Po). Let us construct 
the massless vectors ri and T2 as a linear combination of the two external legs, say Ki and 

K2, 

rf = 1 (i^r + ^K^'^ , r^^ = i (^K^ + ^K^'^ , (2.16) 



with 

P = l-^^1^, and 7 = Ki • K2 + sgn(l, Ki ■ K2)^ {K^ ■ K^f - KfKl .(2.17) 
Then, one builds the niassless vectors t^ and r4 from ri and T2 

^„ _ M^ , r! = (!i!^ , (2^18) 

such that 

rf = , n -rs = Ti •T4 = r2 -Ta = T2 •r4 = n ■ T2 = -T3 ■ u . (2.19) 

The basis Tj can be used to decompose the loop-momentum q, as 

q^ = -p^ + Xirf + X2T^ + XgTg^ + 2;4t|' , (2.20) 

We double this procedure, by introducing a second basis e^ to decompose the other 
loop momentum, 

k^ = -< + yiet + y2e^ + yse^ + y4e^ . (2.21) 

where the basis e^ can be constructed through the same definitions as Tj by replacing 
{Ki, K2) with another couple of external vectors. 

The two basis {tj} and {ej} are adopted for decomposing the solutions of the multiple- 
cuts of the two-loop amplitude. In the following, the decompositions of q and k will be 
uniquely defined by specifying, case by case, po, ti, r2, and tq, ei, 62, respectively. 

2.3 Testing the Integrand-Decomposition 

The integrand-reduction algorithm offers self-consistency checks that ensure the correctness 
of the reconstructed polynomials after the determination of their unknown coefficients. 
These checks, called local and global [N = A^)-tests, are analogous to the tests employed 
in the one- loop integrand-reduction [21,23]. We simply recall here their definitions: 

• The local (N = N)-test monitors the completeness of the polynomial residues, namely 
it can be used to verify that all the necessary ISP's, and their powers, have been prop- 
erly accounted for. Accordingly, after the determination of the polynomial coefficients 
of a given residue, Ajj^...^j^, one can use any solutions of the corresponding ?n,-fold 
cut-conditions other than the ones used for the determination of its coefficients, to 
verify the fulfillment of Eq.( |2.5D . This test can be implemented locally, cut-by-cut, 
and its failure means an incomplete parametrization of the residue. 
We observe that there is only one case where the local {N = N)-test cannot be ap- 
plied: the maximal-cut of a given loop, defined as the cut where all loop momenta are 
frozen by the on-shell cut-conditions. In this case, in fact, the number of solutions of 
the maximal cut is finite, and all of them might have been used for the determination 
of the polynomial coefficients. 

Apart from that, the local {N = N)-test is a very powerful tool for the classification of 
the polynomial structures characterizing the residue of the tti— fold cut for arbitrary 
amplitudes, at any loop, for any m (other than maximal). 



The global {N = N) -test ensures the correctness of the overaU integrand-decomposition. 
Accordingly, after the determination of all polynomial coefficients appearing in the 

must hold for 



r.h.s. of Eq.(2^), the identity between l.h.s. and r.h.s. of Eq.(|27 
arbitrary values of the loop variables. One can therefore verify it: either i) at the 
end of the reduction procedure, or ii) during the reduction, in order to rule out any 
further contribution possibly coming from sub-diagram Mi's. In the latter case, the 
failure of the test indicates that other contributions are missing, and the reduction 
should continue for their detection. 

In the next sections we will present a series of examples that illustrate the integrand- 
reduction algorithm explicitly. The required spinor-algebra has been implemented in Math- 
ematical using the package S@M [45] . 

3. Four-gluon MHV Amplitude in A/" = 4 SYM 

The 4-gluon MHV amplitude in A/" = 4 supersymmetric gauge theory was originally cal- 
culated in Ref. [39]. Two Mi's appearing in the result are the planar and the crossed 
double-box, shown in Fig.Q. 




Figure 1: Two Master Integrals of the 4-gluon MHV amplitude in A/" = 4 SYM: the ladder (left) 
and crossed (right) s-channel double-box. 



3.1 A contribution to the Ladder Amplitude 

Let us consider the contribution to the (leading-color) 4-gluon MHV amplitude in A/" = 4 
SYM coming from the helicity configuration depicted in Figj2|, where only gluons circulate 
in both loops. This case has been discussed in the context of generalised unitarity-based 
method in Ref. [41]. Here we reproduce the same result, and discuss its decomposition in 
terms of the planar double-box MI in Fig.ffl (left), achieved by integrand-reduction. 

On-shell Solutions. The solutions of the 4-point 7fold-cut, Di = . . . = Dq = Dj = 
can be decomposed according to Eqs. ( 2^2^ , 2. 21 ), with the following definitions: 
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In this case, the seven on-shell conditions, Di 
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--P4 



Da = L»7 



A* 

?{7) 












(3.1) 
(3.2) 

0, cannot freeze the 
(3.3) 



10 



4+ 



\ 



/ 



/ 



^'- k 



\ 



Da = (fc-pi)' 
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D7 = (g + fc + p2 



-P3j 



Figure 2: 7fold-cut of the 4-point ladder diagram (s-channel) 

where one of the on-shell cut-conditions imposes a non-hnear relation among X4 and ys, 
which can be imphcitly written as X4 = X4(y3). We choose ya, namely the component of k 
along 63, as the variable parametrizing the infinite set of solutions of the 7fold-cut. 

Residue. The residue of this 4-point 7fold-cut is defined as, 

Al234567(?, k) = Resi234567{^(g, k)} , (3.4) 

where Resi234567{-^(9) k)} is the product of the six 3-point tree-amplitudes sitting in the 
vertices exposed by the cuts. The definition of each 3-point tree-amplitude can be derived 
from the general expression, 

.. (12)=^ A23)V 



/I tree 
^3 



(l-,2-,3^ 



(23)(31)\^(12); 

with a = 0, 1,2 respectively for gluons, fermions, and scalars. In the case at hand, a 
To find the polynomial expression of A1234567, we use the following criteria. 



(3.5) 



0. 



Diagram topology and Vector basis. 

In this 4-point double-box, only three external vectors are independent, therefore 
any possible basis involving them needs to be completed by an additional orthogonal 
vector W3, defined as. 



wS 



{P2 ■ T3)t^ - {P2 ■ T4)t!^ 



{P3 ■ 63)6^ - {P3 ■ (ii)e 



(3.6) 



7"3 • T4 63 • 64 

Irreducible Scalar Products. 

The ISP's which can be formed by the loop variables and the external momenta can 
be chosen to be {q ■ P2) and {k ■ p^). Due to the explicit expression of the on-shell 
solutions, on the cut, one has 

{q ■ P2) -^ {q(j) ■ P2) = (9(7) • ^3) = -3^4(7-4 • P2) , (3.7) 

(fc • ps) -^ {k(j) ■ ps) = (fc(7) • a;3) = ^3(63 • ps) . (3.8) 
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Therefore, A1234567 is as a polynomial in terms of X4 and 2/3, and can be written as 

r r 

Al234567(?, k) = Ci234567,0 + /^ Ci234567,i ("^3 ' if + /_^ Ci234567,i+r- (^3 • kf , (3.9) 

i=l i=l 

where r is the maximum rank in the integration momenta. Due to the orthogonality 
between uj^ and the external momenta, we identify two classes of vanishing integrals. 



4„ / A,^ (u^s-qY 



/"'"/ 



d\ / du 



D, ...D- 



7 



d\ I d'&' ■ ^)" 



0, 
. 



(3.10) 
(3.11) 



Di...Dj 

As a result, the ISP's (a;3 • q) and (073 • k) are spurious, and the polynomial form in Eq.(^Ig|) 
does not generate any additional MI other that the scalar one, whose coefficient is ci234567,o- 

Coefficients. The (2r + 1) unknown coefficients ci234567,j can be determined by sampling 
Eg . (p!^!) and Eq. (|3.9|) on (2r + 1) solutions of the 7fold-cut {q(7),k^Y-^) corresponding to 
(2r + 1) different values of ^3 (remember that 3:4 also depends on 7/3). In this case, we 
finally find that only C1234567 is non- vanishing. 



C1234567,0 = -A*'"'^^(l- , 2", 3+, 4+)s2^S23 , 
Cl234567,i = (1 < i < 2r + 1) . 



(3.12) 
(3.13) 



This means that, after reconstruction, A1234567 is constant, Ai234567(g, ^) = 61234557,0 ; 
which is exactly the result of Ref . [41] . 

3.2 A contribution to the Crossed Amplitude 




D2 = {k+pif 

D3 = {q + k - Pi) 

Ds = (q + paf 

D& = {q-Pif 

Dj = (g + fc + p2 - 



p-i) 



Figure 3: 7fold-cut of the 4-point crossed diagram (s-channel) 

In this example, we consider the contribution to the (subleading-color) 4-gluon MHV 
amplitude in A/" = 4 SYM coming from the helicity configuration depicted in Fig.y, where 
only gluons circulate in both loops. Its decomposition in terms of the crossed double-box 
MI in Fig.|l] (right) is achieved by integrand-reduction. 
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On-Shell Solutions. The solutions of the 7fold-cut, Di = . . . = Dq = Dj = 0, can be 



decomposed according to Eqs.(2.20,2.21), with the fohowing definitions: 

p^ = 0^ e^=P^ e^=P^ (3.14) 

r-o^ = 0^ Tf=p^ T!i=pt. (3.15) 

As before, the on-shell conditions, Di = . . . = Dq = Dj = 0, cannot freeze the loop 
momenta, and one component is left over as a free variable. We choose 2/3, namely the 
component of k along 63 , as the variable parametrizing the infinite set of solutions of this 
7fold-cut. 

Residue. The residue of this 7fold-cut is defined as, 

^1234567iq,k) = ReSi234567{A^(g, fc)} , (3.16) 

where Resi234567{-/V(Q', A;)} is the product of the six 3-point tree-amplitudes sitting in the 
vertices exposed by the cuts, which can be built out of the tree-level expressions in Eg . (p^q) 
used for the planar integrand. 

To find the polynomial expression of A1234567, we use again our two criteria. 

1. Diagram topology and Vector basis. 

As for the planar case, in the considered 4-point double-box, only three external 
vectors are independent, therefore any possible basis involving them needs to be 
completed by an additional orthogonal vector a;3, previously defined in Eq.( |3.6|) . 

2. Irreducible Scalar Products. 

As in the planar case, the ISP's formed by the loop variables and the external mo- 
menta can be chosen to be {q ■ P2) and {k ■ p^). Due to the explicit expression of the 
on-shell solutions, unlike the planar case (see Eq.(p.^)), they both behave linearly in 
1/3. This reflects the fact that these two ISP's are not linearly independent, as one 
can explicitly see using the set of denominators in Fig.^. 

Therefore, A1234567 is a polynomial in 2/3, and can be parametrized as 

r r 

Al234567(5, ^) = X] '=1234567,* (^3 ' ^T = "^ Ci234567,i ((63 ' Ps) Vs)' , (3.17) 

i=0 4=0 

where r is the maximum rank in the integration momenta. The last equation holds because 
(63 • k) = (63 • 64) 2/4 and 7/4 = 0, as required by the seven on-shell conditions. 

Coefficients. One can determine the (r + 1) unknown coefficients 01234567,4, by solving a 
system of equations obtained from sampling Eq.( 3.16|) and Eq.( t3.17|) on (r + 1) solutions 



of the 7fold-cut. The solution of the system finally reads. 



Cl234567,0 


= 


- 


^tree^;^- 


,2 


-,3+,4+) 


^l. 


S22, 


Cl234567,i 


/ 







{i = 


1, 


...,4), 






Cl234567,j 


= 







(5< 


J 


< r) . 







(3.18) 
(3.19) 
(3.20) 
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By using these coefficients in the polynomial expression in Eg. ( 3.17] ) we verify that the 
local-{N = N) test is fulfilled. In practice, we verify the equivalence of Eq.( p.l6 ) and the 
reconstructed polynomial Eq.( p.l7 ) when evaluated in any solution of the 7fold-cut other 
than the ones used to determine the coefficients. 

Notice that in this case the coefficients ci234567,i (i = 1; • • • ;4) are non- vanishing, but 
they do not affect the integrated result, because they multiply spurious integrals that 
vanish upon integration. Nevertheless their presence is important for the completeness of 
the polynomial expression in Eq.( 3.17j ). 



4. The MHV Pentabox in AT = 4 SYM 




D2 = (fc + pa)' 

D3 = {k-pif 

Di = q^ 

D5 = (q + p-if 

De = iq-p4y 

D7 = (q - P4 -ps)^ 

Da = (q + k+p2 + Pa) 



Figure 4: 5-point pentabox diagram. 

In the previous two examples we have discussed the 7fold-cut of a 4-point two-loop 
amplitude. In this section we apply the integrand-reduction to the decomposition in terms 
of Mi's of the 5-point two-loop pentabox in AA = 4 SYM, shown in FigJ^. The expression 
for the pentabox-integrand has been given in a very compact form in Ref. [46]. The 
decomposition involves three types of contributions, coming from: a 5-point 8fold-cut, two 
4-point 7fold-cuts, and two 5-point 7fold-cuts. 

4.1 Five-point Eightfold-Cut 




Figure 5: 5-point 8fold-cut Ai2345678- 
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On-Shell Solutions. The solutions of the 8fold-cut, Di = ... = Dg = 0, depicted in 
Fig.^, can be decomposed according to Eqs.(2.2C,2.21), with the following definitions: 



< = 0^ e^=K, e^=P^ (4.1) 

K = 0^ Tf=p^ T^=p'i. (4.2) 

The eight on-shell conditions Di = ... = Dg = admit four solutions, where the loop 
momenta are completely frozen, namely they are expressed in terms of external kinematic 
variables. 

Residue. The residue of the 8fold-cut is defined as, 

Al2345678(g, k) = ReSl2345678|^^(9, k)j (4.3) 

where the numerator function N can be found in Table I (a) of [46] , and can be parametrized 
as 

^12345678 ('7, k) = Ci2345678,0 + Cl2345678,l {Q ' Pi) + 

+C12345678,2 (k ■ Pi) + Ci2345678,3 {k ' Ps) • (4.4) 

To derive its expression, we considered that this 5-point two-loop integral in four dimensions 
can depend on four external momenta, hence four (out of five) legs can be chosen as a basis 
for decomposing the loop variables. Moreover, this pentabox diagram admits three ISP's, 
and we have chosen {q ■ pi), (k ■ p^), and {k ■ p^). Any other ISP can be expressed as a 
combination of them plus reducible scalar products. 



The definition in Eq.(4.4) is compatible with the existence of four Mi's with denomi- 
nators Di, . . . , Dg: the scalar, plus three rank-1 integrals, each carrying one of the chosen 
ISP in the numerator. 

Coefficients. By sampling N{q, k) at the four solutions of this 8fold-cut we can determine 
the coefficients, and find that ci2345678,o and ci2345678,i are non-vanishing, while ci2345678,2 = 
ci2345678,3 = 0. Therefore, only two (out of four) pentabox-like Mi's will appear in the 
decomposition of the two-loop 5-point amplitude in A/" = 4 SYM. Our result is in agreement 
with the results of Ref. [42,47], although, in order to match the coefficients given in these 
references, one must use an modified expression for A12345678, in which 012345678,1 multiplies 
i{q-p4:-P5) -pi)- 

4.2 Five-point Sevenfold- Cut (i) 

On-Shell Solutions. The solutions of the 7fold-cut, Di = . . . = Dq = Ds = 0, in Fig.| 
can be decomposed according to Eqs.( 2.2C| ,2.21), by using. 



p^ = 0^ e^=P^ e^=P^ (4.5) 

r■o^ = 0^ Tf=p^ T^=P>1. (4.6) 

In this case, the seven on-shell conditions, Di = . . . = Dq = Dg = 0, cannot freeze the 
loop momenta, and the generic solution reads, 

g^t,) = X4T^ , k'^j^ = yse^ , (4.7) 
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1 


) 


n^i = 


...J...' 




i 





■■k 



Figure 6: 5-point 7fold-cut A1234568 

where one of the on-shell cut-conditions imposes a non-hnear relation among X4 and ys, 
which can be imphcitly written as X4 = 234(2/3). We choose 7/3, namely the component of k 
along 63, as the variable parametrizing the infinite set of solutions of the 7fold-cut. 

Residue. The residue of this 7fold-cut is defined as, 

A f h\ T> j N{q,k)-Ai23i5678{q,k)\ 

Al234568W,^) = Resi234568< ^ > ■ (4.8) 

where Ai2345678(9) k) is the polynomial residue of the 8fold-cut, reconstructed in Eq.(^). 
To find the polynomial expression of A1234568) we use the following criteria. 

1. Diagram topology and Vector basis. We observe that this 5-point diagram depends on 
four external momenta and, as in the pentabox case, its integrand does not contain 
any spurious ISP. 

2. Irreducible Scalar Products. The ISP's which can be formed by the loop variables 
and the external momenta can be chosen to be {q ■ pi), (q ■ P2), [k ■ p^), {k ■ p^). Due 
to the explicit expression of the on-shell solutions, on the cut, one has 

{q ■ pi) -^ (g(7) • pi) = -a;4(r4 • pi) , (4.9) 

{q ■ P2) -^ {q(j) ■ P2) = -Xi{Ti ■ P2) , (4.10) 

{k ■ P3) -^ (fc(7) • ps) = 2/3(63 • ps) , (4.11) 

{k ■ P4) -^ (/c(7) • pa) = 2/3(63 • Pa) ■ (4.12) 

Therefore, A1234567 is polynomial in X4 and 2/3, and can be written as 

r r 

Ai234568(Q, k) = Ci234568,0 + ^ Ci234568,j {P2 'QY + '^ Ci234568,i+r (^53 ' k)' , (4.13) 

i=l i=l 

where r is the maximum rank in the integration momenta. 

Coefficients. As by-now understood, the unknowns 01234568,4 are found by sampling 
Eg. ([4.^) and Eg. (4. 13]) on (2r -|- 1) solutions of the 7fold-cut. Also in this case we find 



that Ai234568(Q) ^) is a trivial polynomial, namely just a constant, because only 01234568,0 
is non- vanishing. 

Let us finally remark that integrands with numerator (p2 • qY, and (p3 • /c)* are non- 
vanishing, and would be Mi's. They do not show up in this case because they are multiplied 
by null coefficients. 
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4.3 Five-point Sevenfold-Cut (ii) 



pA = 


...J.. J 




J 




A 







'■k 



Figure 7: 5-point 7fold-cut A1234678 

This case is specular to the previous one and can be treated by relabeling the external 
momenta (1 -f-)- 2, 3 -f-)- 5). As before only ci234678,o is non-vanishing, therefore only the 
scalar integral contribute. 

4.4 Four-point Sevenfold-Cut (i) 



\ 








i 






i 





Figure 8: 4-point 7fold-cut A1234578 

The 7fold-cut in Fig. ^ is equivalent to the ladder case treated in Sec. |3.1| . The solutions 
of Di = . . . = D^ = Dt = Ds = can be decomposed according to Eqs.( p.2d , 2.21 ), by 
using 






0>' , 
0^ , 



Pi 
= ^3 



P2 ' 

- P^ _ '^45 ^/j 

" -^45 9 p ^ ^1 
2i-45 • Ti 



(4.14) 
(4.15) 



As in Sec. 3.1, the seven on-shell conditions are not sufficient to freeze the loop momenta. 



and one component is left over as a free variable. We choose 2/3, namely the component of 
k along 63, as the variable parametrizing the infinite set of solutions. 

The residue of the 7f old-cut, Di = . . . = D^ = Dj = Dg = 0, is defined as. 



Al234578(g, ^) = Res 



1234578 ' 



N{q,k)-A 



12345678 



{q,k) 



De 



(4.16) 



where Ai2345678('7) ^) is the polynomial residue of the 8fold-cut, reconstructed in Eq.(^). 
The polynomial expression of A 1234578 is equivalent to the one given in Eq.(^^), where Wg , 



defined in Eq.(3^), must be constructed with the basis vector {e,} used in Eq.( [4.15| ). 

We determine the unknown coefficients and find that only 01234578,0 is non-vanishing. 
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g-p4-t---:-- + ---i---T-fc 

/ 

4 2 

3 

Figure 9: 4-point 7fold-cut A1235678 



4.5 Four-point Sevenfold-Cut [ii) 



Also the case in Fig.^ falls in the two-loop 4-point category discussed in Sec. 3.1 and in 



the previous section. The solutions of the 7fold-cut, Di = . . . = D3 = Z^s = . . . = Dg = 0, 
can be decomposed according to Eqs. (2^2^,2. 21), with the following definitions: 



P^ = 0^ e^=P^ e^=P^ (4.17) 

-o" = -pr -r=p^ < = pt.-^^^r{- (4.i8) 

Again, we find it convenient to use 2/3 as free variable to parametrize the infinite set of 
solutions. 

The residue, defined as 

A / /\ T3 f A^(9,^)- Ai2345678(9,fc)l . 

Al235678W,A;) = Resi235678< ^ > , (4.19) 



with Ai2345678(9) ^) being the the 8fold-cut polynomial in Eq.( [4.4[ ), can be written as in 
Eg. (^^) , where Wg must be replaced by 

. = (P5 ■ e3)e^ - (P5 • e,)el .4 20) 

^ 63-64 

By polynomial sampling we find that only 61234578,0 is non-vanishing. 

4.6 Reconstructed Integrand 

Combining the results of the previous sections, the numerator of the planar 5-point pentabox 
diagram can be decomposed as, 

lSS{q,k) = Ai2345678(5,^) + 

+Ai234568(9, k)D7 + Ai234578(9, k)De + 
+Ai234678(g, k)D5 + Ai235678(9, k)Di = 
= Ci2345678,0 + Cl2345678,l (Q ' Pi) + 
+£1234568,0-^7 + £1234578,0-^*6 + 
+£1234678,0-^5 + £1235678,0-^4 , (4-21) 

which corresponds to a decomposition of the integral in terms of six Mi's, 
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Cl2345678,0 4 



+Cl234568,0 




+C1234678,0 



+ Ci2345678,l 4 



1 4 

N 

+ Ci234578,0 



2 3 

1 5 

N 

+ Cl235678,0 




+ 



+ 



(4.22) 



We checked the correctness of the Eq.( [4.21 ) through the global-{N = N) test, namely by 
verifying the identity of the l.h.s and of the r.h.s. for arbitrary values of q and k. 



5. The MHV Pentacross in A/" = 4 SYM 




D2 = (k + p2f 

D3 = (fc + (7-p4 -Pif 

D4 = q^ 

D5 = (q + Psf 

De = iq-p4f 

D7 = (g -P4 -ps)^ 

-Ds = {q + k+p2 +P3)^ 



Figure 10: 5-point Pentacross (non-planar) 

In this section we present the result of the integrand reduction of the MHV pentacross 
diagram in J\f = 4 SYM, depicted in Fig. 10. The expression for the its integrand has been 
given in Table I (b) of Ref. [46], and happens to have the same expression of the planar 
diagram, although it is sitting on a different set of denominators, listed in FigJTo|. 

The integrand-reduction follows the same pattern as described in the Section 3 and 4. 
In particular the expressions of the polynomial residues of the 5-point 8— fold cut, and 4- 
point 7— fold cut, can be obtained following the same procedures as in Sec. iA and Sec. |3.2| , 
respectively. The 5-point 7— fold cut deserves a dedicated discussion. 
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5.1 Five-point Sevenfold-Cut 

On-Shell Solutions. The solutions of the 5point 7fold-cut, Di = . . . = Dq = Dg = 0, 



depicted in Fig. 11 can be decomposed according to Eqs.( 2.20 , 2.21 ), with the following 
definitions: 

(5.1) 
(5.2) 

The on-shell conditions, Z?i = . . . = Dq = Dg, = 0, cannot freeze the loop momenta, and 
we choose 2/3 to parametrize the infinite set of solutions. 



K = o^ 


Gi — Pi , 


e^=P^ 


r-(; = o^ 




<=P^ 




3 
Figure 11: 5-point 7fold-cut A1234568 

Residue. The residue is defined as, 

A f l\ T) j N{q,k)- Ai23i5678iq,k)\ , . 

Al234568(g,«:) = Resi234568< ^ > ■ (5.3) 

The function A12345678 is the 8fold-cut polynomial of the pentacross. To find the polynomial 
expression of A1234568) we use the following criteria. 

1. Diagram topology and Vector basis. We observe that this 5-point diagram depends 
on four external momenta, therefore its integrand does not contain any spurious ISP. 

2. Irreducible Scalar Products. The ISP's formed by the loop variables and the external 
momenta can be chosen to be (q ■ pi), {q ■ P2), {k ■ ps), {k ■ p^). Due to the explicit 
expression of the on-shell solutions, they all behave linearly in y^. This reflects the 
fact that although irreducible, these four ISP's are not linearly independent. 

Therefore, A1234568 is as a polynomial in 2/3, and can be written as 

r 

Al234568(g, k) = Ci234568,0 + ^ Ci234568,j ipi ' qY , (5-4) 

1=1 

where r is the maximum rank in the integration momenta. 

Coefficients. As by-now understood, the unknowns 01234568,2 are found by sampling 
Eq.(^.3|) and Eq.(^.4|) on (r + 1) solutions of the 7f old-cut. Also in this case we find 
that Ci234568,0 is the only non-vanishing coefficient. 

Let us finally remark that integrands with numerator (pi • qY would generate Mi's. 
They do not show up in the final result because they are multiplied by null coefficients. 
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5.2 Reconstructed Integrand 

The integrand-decomposition for the MHV pentacross diagram given in Table I (b) of 
Ref. [46] has the same structure of the planar pentabox, 

N{q,k) = Ai2345678(g, fc) + 

+Ai234568(9, k)D7 + Ai234578(9, k)De + 
+Ai234678(9, k)D5 + Ai235678(9, k)Di = 
= Ci2345678,0 + Cl2345678,l {Q ' Pi) + 
+£1234568,0-^7 + Ci234578,0-D6 + 
+£1234678,0-^5 + 01235678,0-^4 , (5.5) 



and contains the same coefficients of Eq.( 4.2l| ). The above decomposition has been verified 
to fulfill the glohal-{N = N) test. Therefore, the final expression of the pentacross diagram 
in terms of Mi's reads, 




Cl2345678,0 4 



+Cl234568,0 




+ Ci2345678,l 4 V(<?-pi) 



1 4 

+ Ci234578,0 




+ 



(5.6) 



+C1234678,0 



+ Cl235678,0 




6. Conclusions 

We illustrated a first implementation of the integrand-reduction method for two-loop scat- 
tering amplitudes. We have shown that the residues of the amplitudes on the multi-particle 
cuts are polynomials written in terms of independent irreducible scalar products formed by 
loop momenta and either external momenta or polarization vectors built out of them. The 
independence conditions among irreducible scalar products can be investigated through 
their polynomial behavior in terms of the components of the loop momenta still undeter- 
mined after imposing the on-shell cut-conditions. 

The reduction of the amplitudes in terms of master integrals can be realized through 
polynomial fitting of the integrand, without any need of an apriori knowledge of the integral 
basis. We discussed how the polynomial shapes of the residues determine the basis of master 
integrals appearing in the final result. In particular, we have found that the multiparticle 
residues of amplitudes with less then five external legs can eventually be written in terms 
of spurious irreducible scalar products which do not generate any master integral upon 
integration. 
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We applied the integrand-reduction algorithm to cases of modest complexity, such as 
planar and non-planar contributions to the 4-point MHV and 5-point MHV amplitudes in 
N = A SYM. We worked out the polynomials parametrizing the residues at the 4-point 
7fold-cut, 5-point 7fold-cut, 5-point 8fold-cut which contribute to the decomposition of the 
considered amplitudes. 

In this work we identified general principles which could guide in the classification 
of all polynomial residues required by the complete integrand-reduction of arbitrary two- 
loop amplitudes. The technique we have presented extends the well-established analogous 
method for one-loop amplitudes, and can be considered a preliminary study towards the 
systematic reduction at the integrand-level of multi-loop amplitudes in any gauge theory, 
suitable for their automated semianalytic evaluation. 
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